Measurements of higher order noise correlations in a quantum dot 
with a finite bandwidth detector 
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We present measurements of the fourth and fifth cumulants of the distribution of transmitted 
charge in a tunable quantum dot. We investigate how the measured statistics is influenced by the 
finite bandwidth of the detector and by the finite measurement time. By including the detector 
when modeling the system, we use the theory of full counting statistics to calculate the noise levels 
for the combined system. The predictions of the finite-bandwidth model are in good agreement with 
measured data. 
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Current fluctuations in mesoscopic systems have been 
extensively studied due to the extra information they give 
in comparison to measurements of the mean current 
The focus has traditionally been on investigations of the 
shot-noise, which for classical systems arises due to the 
discreteness of the electron charge. The theory of full 
counting statistics (FCS) was introduced as a new way 
of examining current fluctuations Q. With the FCS, fluc- 
tuations are studied by counting the number of electrons 
that pass through a conductor within a fixed period of 
time. This gives direct access to the probability distri- 
bution function p to (N), which is the probability that TV 
electrons are transferred within a time interval of length 
to. From the distribution function, not only the shot 
noise but also correlations of higher order can be calcu- 
lated. 

The third moment of a tunneling current has been 
shown to be independent of the thermal noise [!, 0|, 
thus making it a potential tool for investigating electron- 
electron interactions even at elevated temperatures. 
Higher order moments in strongly interacting systems 
are predicted to depend heavily on both the conductance 
Q and on the internal level structure [|| of the system. 
Determining higher order moments may therefore give a 
more complete characterization of the electron transport 
process. This can be of importance for realizing measure- 
ments of electron correlation and entanglement effects in 
quantum dots @, HI • 

Experimentally, the third moment of the current distri- 
bution function has been measured for a tunnel junction 
[§| as well as for a single quantum dot (QD) [l(| E| and 
a double QD Q. In quantum optics, higher order mo- 
ments are routinely measured in order to study entangle- 
ment and coherence effects of the electromagnetic field 
Here, we set out to measure the fourth and fifth cu- 
mulant of the distribution function for charge transport 
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through a QD. 

In general, experimental measurements of FCS for elec- 
trons are difficult to achieve due to the need of a sensi- 
tive, high-bandwidth detector capable of resolving indi- 
vidual electrons [HI, 0, [HI]. However, a more funda- 
mental complication with the measurements is that most 
forms of the FCS theory assume the existence of (1) a 
detector with infinite bandwidth and (2) infinitely long 
data traces. Since no physical detector can fulfill these 
requirements, every experimental realization of the FCS 
will measure a distribution which is influenced by the 
properties of the detector. Here, we investigate how the 
violation of the two assumptions modifies the measured 
statistics. By including the detector in the model, the 
FCS for the combined QD-detector system can be calcu- 
lated. This model can explain the results for higher order 
cumulants measured with a finite bandwidth detector. 

The sample consists of a QD [dotted circle in the in- 
set of Fig. [ljb)] with a nearby quantum point contact 
(QPC) used for reading out the charge state of the QD 
[1 61 ] . The structure was fabricated using scanning probe 
lithography [13] on a GaAs/Alo.3Gao.7As heterostructure 
with a two-dimensional electron gas 34 nm below the sur- 
face. The gates G\ and G2 are used to tune the height of 
the tunneling barriers connecting the dot to source and 
drain leads, while the P-gate is used to keep the conduc- 
tance through the QPC in a regime where the sensitiv- 
ity to changes in its electrostatic potential is maximal. 
All measurements were performed in a dilution refriger- 
ator with a base temperature of 20 mK. The electronic 
temperature extracted from the width of Coulomb block- 
ade resonances measured in the low bias regime [l8j was 
190 mK. 

Due to the electrostatic coupling between the QD and 
the QPC, the addition of an extra electron on the dot will 
cause a change in the QPC conductance. By performing 
time-resolved measurements of the current through the 
QPC, tunneling of single electrons can be detected in 
real-time 

0, [Hill- In this experiment, the QPC was 
voltage biased with Vqpc = 250 /iV '. The current signal 
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was sampled at 100 kHz, software filtered at 4 kHz using 
a 8th order Butterworth filter and finally resampled at 
20 kHz to keep the amount of data manageable. 

To measure the higher cumulants for the current 
through the QD, one has to generate the experimental 
probability density function p ta (N). This is done by 
splitting a time trace of length T into m — T/to in- 
tervals and counting the number of electrons entering 
the dot within each interval. The cumulants are then 
calculated directly from the distribution function. In a 
previous experiment on a similar system the second and 
third cumulants were measured [lOj . To extend the anal- 
ysis to higher cumulants, it is necessary to increase the 
length of the time traces in order to collect more statis- 
tics. Here, we present cumulants extracted from time 
traces of length T = 10 minutes. The quality of the 
data allows to measure up to the fifth cumulant, which is 
two orders higher than reported in previous experiments 
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FIG. 1: (a-d) Normalized cumulants C n /Ci versus dot asym- 
metry, a = (Ti n — r out ) /(Ti n + r out ). The solid lines are theo- 
retical predictions assuming a perfect detector, C2/C1 — (1 + 
a 2 )/2, C3/C1 = (l + 3a 4 )/4, C4/C1 = (1 + a 2 - 9a 4 + 15a 6 )/8 
and C5/C1 = (1 + 30a 4 - 120a 6 + 105a 8 )/16. The dashed 
lines show the cumulants calculated from the model defined 
by Eq. fTJ in the text. The inset in (b) shows the quantum 
dot with integrated charge read-out used in the experiment. 
The inset in (c) shows the variation of the total tunneling rate 
Ttot — Tin + Tout for the different measurement points. 

The results are shown in Fig. [TJ where we plot the nor- 
malized cumulants for different values of the asymmetry 
of the tunneling rates, a = (T- m — r ou t)/(r in +r out ). Here, 
Tin and r out are the rates for tunneling into and out of 
the dot, respectively. The asymmetry is tuned by shifting 
the voltage on gate G\ by an amount AV and at the same 
time applying a compensating voltage — AV on gate G2. 
With the two gates having a similar lever arm on the dot, 
the electrochemical potential of the QD remains at the 
same level, but the height of the tunneling barriers be- 
tween the dot and the source and drain leads will change. 



Doing so, we could tune the asymmetry from a = —0.94 
to a = +0.25 while still keeping both tunneling rates 
within the measurement bandwidth and avoiding charge 
rearrangements. To get data for the full range of asym- 
metry, we did a second measurement at a different gate 
voltage configuration. For the second set of data, the 
asymmetry was tuned from a = 0.07 to a — 0.93. The 
stars and the circles in Fig. [T] represent data from the 
two different sets of measurements. The measurements 
were performed with a QD bias of Vbi as = 2.5 mV, with 
the electrochemical potential of the dot far away from 
the Fermi levels of the source and drain leads. This is 
to ensure that tunneling due to thermal fluctuations is 
sufficiently suppressed 

The solid lines in Fig. [1] depict the theoretical predic- 
tions calculated from a two-state model [2l[ . The analyti- 
cal expressions are given in the figure caption. The higher 
cumulants show a complex behavior as a function of the 
asymmetry, with local minima at a — ±0.6 for C4/C1 
and at a = ±0.8 for C5/C1. The fifth cumulant even 
becomes negative for some configurations. The experi- 
mental data qualitatively agrees with the theory, but for 
small values of the asymmetry there are deviations from 
the expected behavior. The deviations are stronger for 
the first set of data (stars). Since the tunneling rates in 
the first measurement was about a factor of three higher 
than in the second measurement [see inset of Fig. [He)], 
we suspect the finite bandwidth of the detector to be a 
possible reason for the discrepancies. 

Recently, Naaman et al. 22] pointed out that mea- 
surements of the transition rates of a Poisson two-state 
system using a finite bandwidth detector always leads 
to an underestimate of the rates. To determine the 
rates correctly, the detection rate Tdot of the detector 
must be taken into account. In the low-bias, weak cou- 
pling Coulomb blockade regime, the QD can be modeled 
as a Poisson two-state system. The two states corre- 
spond to zero or one excess electron on the dot, and 
the transitions between the two states occur whenever 
an electron tunnels into or out of the dot. The proba- 
bility distribution for the times needed for an electron 
to tunnel into or out of the QD follows the exponential 

Pin/out (*) = r in/out exp(-r in/out t) [m. 

An example of a probability distribution taken from 
measured data is shown in Fig. [Ha). The long-time be- 
havior is exponential, but for times t < 100 fis there is a 
sharp decrease in the number of counts registered by the 
detector. From the figure, we can estimate Tdct, which 
is the average time it takes for the detector to register 
an event. We find Td c t = 70 /is, giving a detection rate 
of Tdct = l/^dct = 14 kHz. Note that the detection rate 
Tdct does not only depend on the measurement band- 
width but also on the signal-to-noise ratio of the detec- 
tor signal as well as the redundancy needed to minimize 
the risk of detecting false events [23| . All tunneling rates 
presented in the following have been extracted from dis- 
tributions such as the one shown in Fig. [D(a), using the 
methods described in Ref. [13] with Tdct = 14 kHz. 
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FIG. 2: (a) Probability density of time needed for an electron 
to tunnel into the dot. Note the sharp decrease in counts for 
t < 100 /j,s due to the finite bandwidth of the detector. The 
black curve is a long-time exponential fit with F = 1.39 kHz. 
(b) Model for the dot-detector system. A state (n, m) cor- 
responds to n electrons on the dot while the detector at the 
same time is measuring m electrons, (c) Higher cumulants 
versus relative detection bandwidth Tdct/tTin + r out ), calcu- 
lated from the model in (b). The cumulants are normalized 
to the results from the infinite-bandwidth case. The influ- 
ence of the finite bandwidth is maximal when the asymmetry 
a = (Ti n - r ut)/(r in + r out ) is zero. 



The finite bandwidth will also influence the FCS mea- 
sured by the detector. Following the ideas of Ref. [12], we 
account for the finite bandwidth by including the states 
of the detector in the model. Figure HJb) shows the four 
possible states of the combined dot-detector model. The 
state (n + 1, n) refers to a situation where there are n + 1 
electrons on the dot, while the detector at the same time 
reads n electrons. The transition from the state (n+1, n) 
to the state (n + l,n+ 1) occurs when the detector reg- 
isters the electron. This process occurs with the rate of 
the detector, Tdet- 

To calculate the FCS for the QD-detector system, 
we write the master equation P = MP, with P = 
[(n, n), (n + 1, n), (n, n + 1), (n + 1, n + 1)] and 
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In the above matrix, we have included the counting fac- 
tor e lx at the element where the detector registers an 
electron tunneling into the dot [see dashed circle in Fig. 
[2Kb)] . The statistics obtained in this way relates di- 
rectly to what is measured in the experiment. Using 
the methods of Ref. [2l|, we calculate the first few cu- 
mulants for the above expression as a function of rela- 
tive bandwidth k — rd t/(ri n + r out ) and asymmetry 



a = (l^in — r ou t)/(rin + r out ). The normalized second 
and third cumulants take the form 
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In Fig. [2Kc) we plot the second and third cumulants 
from Eq. and Eq. for different values of asymme- 
try a and relative bandwidth k. The cumulants have been 
normalized to the values for the infinite bandwidth de- 
tector. With Tdot 3> r in + r out , the cumulants approach 
the infinite bandwidth result, as expected. However, even 
with Tdct = 10(rin + r out ) and perfect symmetry (a = 0), 
the second cumulant deviates by almost 10% and the 
third cumulant by more than 20% from the perfect de- 
tector values. As the bandwidth is further decreased, 
the deviations grow stronger and reach a maximum as 

Tdct = Tin H~ r ou t . With 

Tdot Tin + Tout, the cumulants 
once again approach the perfect detector values. When 
the detector is much slower than the underlying tunnel- 
ing process, it will only sample the average population of 
the two states. In this limit, the dynamics of the system 
does not interfere with the dynamics of the detector and 
we recover the correct relative noise levels. It should be 
noted that this is true only for the noise relative to the 
detected mean current. Since the detector will miss most 
of the tunneling events, the absolute values of both the 
current and the noise will be underestimated. 

We have also performed the analysis for the fourth and 
fifth cumulants. We do not show the analytical expres- 
sions due to space limitations; however, the results cor- 
responding to the experimental configuration are repre- 
sented by the dashed lines in Fig. 1. Over the full range 
of bandwidth and asymmetry, we find that the noise de- 
tected with the finite bandwidth system is always lower 
than for the ideal detector case. The reduction can be 
qualitatively understood by considering the probability 
distribution p to (N). The finite bandwidth makes it less 
probable to detect fast events, meaning that the prob- 
ability of detecting a large number of electrons within 
the interval to will decrease more than the probability 
of detecting few electrons. This will cut the high-count 
tail of the distribution and thereby reduce its width (C2) 
and its skewness (C3). An interesting feature is that the 
cumulants calculated for a less symmetric configurations 
[a = 0.9 in Fig. HJc)] show less influence of the finite 
bandwidth. 

A second limitation of a general FCS measurement is 
the finite length of each time trace. In order to generate 
the experimental probability density function p to (iV), the 
total trace of length T must be split into m = T/tg inter- 
vals, each of length to. Most FCS theories only predict 
results for the case to S> 1/r, where T is a typical transi- 
tion rate of the system. In the experiment, it is favorable 



4 



to make to as short as possible in order to increase the 
number of samples m = T/t . This will improve the 
quality of the distribution and help to minimize statisti- 
cal errors. 

The condition to 3> T is imposed by the approximation 
that the cumulant generating function (CGF) S(x) f° r 
Pt (N) only depends on the lowest eigenvalue A m ; n of 
the master equation matrix M x , with S{\) = — toA m i n . 
A FCS valid for finit e to must include all eigenvalues and 
eigenvectors of M x [21| • The corresponding expression 
is 

exp[S( X )] = (qo\ P [n) ) exp(-t A n )(gW|po), (4) 

where {q^\ and \p^) are the left and right eigenvectors 
of the matrix M x , A n are the eigenvalues of M x and 
(<lo |; bo) are the eigenvectors corresponding to the lowest 
eigenvalue A m i n . The cumulants generated from the CGF 
in Eq. ([4]) will in general be a function of to- 




FIG. 3: Normalized cumulants evaluated for different lengths 
of the time interval to. The symbols show the experimental 
data, extracted from a time trace of length T = 10 minutes, 
containing 350595 events, with a = 0.053, and r to t = 3062 Hz. 
The solid lines are calculations from the FCS given by Eq. 
@ in the text, while the dashed lines are the asymptotes 
for to — > tx). The inset shows a magnification of the vertical 
axis (horizontal axis unchanged) for C4/C1 and C5/C1 for 
(N) > 0.6. 

To investigate how small to can be before systematic 
errors become relevant, we calculate the cumulants from 
the CGF of Eq. (j4|) with the master equation matrix M x 
of Eq. JTJ). The results are shown in Fig. [31 where we 
plot the normalized cumulants as a function of the mean 
number of counts per interval, (N) = to/(l/ri n + l/r out ). 
The symbols show cumulants extracted from measured 
data (T = 10 minutes, a = 0.053, T in + r out = 3062 Hz 
and Tdct = 14 kHz), while the solid lines are results from 
the CGF for the same set of parameters. The dashed 
lines are the asymptotes for the limiting case to 00. 



In general, the data and the theory are in good agree- 
ment. There are some deviations in the fourth and 
fifth cumulants for large to ((N) > 6 in Fig. [3]), but 
these are statistical errors due to the shortness of the to- 
tal time trace. For short to, all cumulants converge to 
C n /Ci — » 1. This is because as (N) <C 1, the probability 
distribution p to (N) will be non-zero only for N = and 
N = 1, withp to (0) = 1-q, p t0 (l) = q and q = (N). This 
is the definition of a Bernoulli distribution, for which the 
normalized cumulants C n jC\ — > 1 as q — > [HI • 

Focusing on the other regime, (N) > 1, we see that cu- 
mulants of different orders converge to their asymptotic 
limits for different values of to- The second cumulant 
needs a longer interval to to reach a specified tolerance 
compared to the higher cumulants. This is of interest 
for the experimental determination of higher cumulants. 
By choosing a shorter value of to when calculating higher 
cumulants, the amount of samples m = T '/to can be in- 
creased. For the data in Fig. [TJ the cumulants were cal- 
culated with intervals to giving (N) = 15 for C2, (N) = 6 
for C 3 , (N) = 3 for C 4 and (N) = 2 for C 5 . The maximal 
deviations between the correct cumulants and the ones 
determined with a finite length to can be estimated by 
checking the convergence for all values of the asymmetry. 
For the data shown in Fig. [2 we find AC2/C1 = 0.007, 
AC3/C1 = 0.009, AC4/C1 = 0.01 and AC5/C1 = 0.03. 

Coming back to the results of Fig. [TJ we are now able 
explain why the measured cumulants show lower values 
compared to the perfect-detector theory. The dashed 
lines in Fig. [JJare the cumulants calculated from the com- 
bined QD-detector model of Eq. ([T]), with T^t = 14 kHz. 
The overall agreement is good, especially since no fitting 
parameters are involved. Higher cumulants end up to be 
slightly lower than theory predicts. We speculate that the 
deviations could be due to low-frequency fluctuations of 
the tunneling rates over the time of measurement. 

In conclusion, we have measured the first five cumu- 
lants of the distribution of charge transmitted through a 
QD. The ability to measure higher cumulants shows that 
we can determine the distribution function very precisely. 
The high accuracy of the technique makes it a promis- 
ing tool for probing subtle effects in the transport statis- 
tics of more complex QD systems. We have found that 
the measured statistics depends strongly on the band- 
width of the charge detector. By including the detector 
in the model, we show that the framework of FCS can 
be used to predict noise levels for systems with a finite 
bandwidth detector. The principle is general and can be 
applied to any rate-equation model used for calculating 
the FCS. Financial support by the NCCR Nanoscience 
through the Swiss Science Foundation (Schweizerischer 
Nationalfonds) is gratefully acknowledged. 



■5 



[1] Y. M. Blanter and M. Biittiker, Physics Reports 336, 1 
(2000). 

[2] L. S. Levitov, H. W. Lee, and G. B. Lesovik, J. Math. 

Phys. 37, 4845 (1996). 
[3] T. Fujisawa, T. Hayashi, R. Tomita, and Y. Hirayama, 

Science 312, 1634 (2006). 
[4] L. S. Levitov and M. Reznikov, Phys. Rev. B 70, 115305 

(2004). 

[5] D. Bagrets, Y. Utsumi, D. Golubev, and G. Schon, 
Fortschritte der Physik 54, 917-938 (2006). 54, 917 
(2006). 

[6] W. Belzig, Phys. Rev. B 71, 161301(R) (2005). 
[7] D. Loss and E. V. Sukhorukov, Phys. Rev. Lett. 84, 1035 
(2000). 

[8] D. S. Saraga and D. Loss, Phys. Rev. Lett. 90, 166803 
(2003). 

[9] B. Reulet, J. Senzier, and D. E. Prober, Phys. Rev. Lett. 
91, 196601 (2003); Y. Bomze et al, Phys. Rev. Lett. 95, 
176601 (2005). 

[10] S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, 

T. Ihn, P. Studerus, K. Ensslin, D. C. Driscoll, and A. C. 

Gossard, Phys. Rev. Lett. 96, 076605 (2006). 
[11] S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, 

P. Studerus, T. Ihn, K. Ensslin, D. C. Driscoll, and A. C. 

Gossard, Phys. Rev. B 74, 195305 (2006). 
[12] L. Mandel and E. Wolf, Optical Coherence and Quantum 

Optics (Cambridge University Press, 1995). 
[13] W. Lu, Z. Ji, L. Pfeiffer, K. W. West, and A. J. Rimberg, 

Nature 423, 422 (2003). 
[14] T. Fujisawa, T. Hayashi, Y. Hirayama, H. D. Cheong, 



and Y. H. Jeong, Appl. Phys. Lett. 84, 2343 (2004). 
[15] J. Bylander, T. Duty, and P. Delsing, Nature 434, 361 
(2005). 

[16] M. Field, C. G. Smith, M. Pepper, D. A. Ritchie, J. E. F. 

Frost, G. A. C. Jones, and D. G. Hasko, Phys. Rev. Lett. 

70, 1311 (1993). 
[17] A. Fuhrer, A. Dorn, S. Liischer, T. Heinzel, K. Ensslin, 

W. Wegscheider, and M. Bichler, Superl. and Microstruc. 

31, 19 (2002). 

[18] L. P. Kouwenhoven, C. M. Marcus, P. M. McEuen, 
S. Tarucha, R. M. Westervelt, and N. S. Wingreen, in 
Mesoscopic Electron Transport, edited by L. L. Sohn, 
L. P. Kouwenhoven, and G. Schon (Kluwer, Dordrecht, 
1997), NATO ASI Ser. E 345, pp. 105-214. 

[19] R. Schleser, E. Ruh, T. Ihn, K. Ensslin, D. C. Driscoll, 
and A. C. Gossard, Appl. Phys. Lett. 85, 2005 (2004). 

[20] L. M. K. Vandersypen, J. M. Elzerman, R. N. Schouten, 
L. H. Willems van Beveren, R. Hanson, and L. P. 
Kouwenhoven, Appl. Phys. Lett. 85, 4394 (2004). 

[21] D. A. Bagrets and Y. V. Nazarov, Phys. Rev. B 67, 
085316 (2003). 

[22] O. Naaman and J. Aumentado, Phys. Rev. Lett. 96, 
100201 (2006). 

[23] C. E. Shannon, Proc. Institute of Radio Engineers 37, 10 
(1949). 

[24] E. W. Weisstein, "Bernoulli Distribution." 
From MathWorld - A Wolfram Web Resource. 
http://mathworld.wolfram.com/BernoulliDistribution.html. 



